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Abstract 

A new improved quark mass density- dependent model including u, d quarks, a mesons, lo mesons 
and p mesons is presented. Employing this model, the properties of nuclear matter, neutron matter 
and neutron star are studied. We find that it can describe above properties successfully. The results 
given by the new improved quark mass density- dependent model and by the quark meson coupling 
model are compared. 
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I. INTRODUCTION 



In our previous papers [1-6], a new quark meson coupling model bases on quark mass 
density- dependent(QMDD) model is presented. The QMDD model suggested by Fowler, 
Raha and Weiner [7] firstly assuming that the masses of u, d and s quarks (and the corre- 
sponding antiquarks) satisfy: 

B - 
mq^- — {i ^ u,d,u,d) (1) 

B 

ms,s = rriso + - — (2) 
3ns 

where is the baryon number density, m^o is the current mass of the strange quark, and 
B is the bag constant. As was explained in Refs.[l, 2, 5, 8], the ansatz Eqs. (1) and 
(2) corresponds to a quark confinement hypothesis because when y — > oo, — > and 
ruq ^ oo, it prevents the quark goes to infinity. It is shown that the properties of strange 
quark matter in the QMDD model are nearly the same as those obtained in the MIT bag 
model [9, 10]. But the basic difference is that instead of the MIT bag boundary condition, 
we have the density- dependent masses of quarks in QMDD model according to Eqs. (1) 
and (2). It means that the ansatz Eqs. (1) and (2) can replace MIT bag boundary condition 
and get the nearly the same results. 

Quark- meson coupling(QMC) model suggested by Guichon [11] firstly is a famous hybrid 
quark meson model which can describe many physical properties of nuclear matter and nuclei 
successfully [12]. In this model, the nuclear system was suggested as a collection of MIT 
bag and mesons. The interactions between quarks and mesons are limited within the MIT 
bag regions. As was pointed in Refs. [1, 2, 6], this model has two major shortcomings: 
(1) It is a permanent quark confinement model because the MIT bag boundary condition 
cannot be destroyed by temperature and density. Therefore, it cannot describe the quark 
dcconfincmcnt phase transition. (2) It is difficult to do nuclear many-body calculation 
beyond mean field approximation(MFA) by means of QMC model, because wc cannot find 
the free propagators of quarks and mesons easily. The reason is that the interactions between 
quarks and mesons are limited within the bag regions, the multirefiection of quarks and 
mesons by MIT bag boundary must be taken into account for getting the free propagators. 
These two shortcomings come from MIT bag constrain all. 
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To overcome these two shortcomings, we suggested an improved quark mass density- de- 
pendent(QMDD) model in Refs. [1-6]. We added the a- meson and u- meson to improve 
the QMDD model. Instead of the MIT bag, after introducing the nonlinear interaction of 
(j-mesons and qqc coupling, we construct a Friedberg- Lee soliton bag in nuclear system. 
The quark masses are still density- dependent. The interactions between quarks and mesons 
are extended to the whole system. Since the MIT bag constraint is given up, our improved 
QMDD(IQMDD) model can describe the quark deconfinement phase transition [6] and do 
the nuclear many- body calculations beyond MFA in principle. We have proved that our 
model can successfully describe the saturation properties, the equation of state, the com- 
pressibility and the effective nucleon mass of symmetric nuclear matter and give a reasonable 
critical temperature of quark deconfinement [1-6]. 

The motivation of this paper is to extend our study to asymmetric nuclear matter, es- 
pecially to the neutron matter and the neutron star. It means that we must consider the 
isosipn dependence and distinguish the u-quark and d-quark. We will add the isospin vector 
p mesons to improve the IQMDD model in this paper. We hope to compare the results of 
IQMDD model with those obtained by QMC model and QHD-II model for neutron matter 
and neutron star. In order to find their differences and similarity explicitly, we will use the 
same approximation as that of the QMC model [13] in our calculations. Though the study 
of neutron star employing QMC model in Ref. [14] is too simple, but in order to exhibit 
the basic differences between the IQMDD model and the QMC model, we still consider the 
neutron star by using the same approximation as Ref. [14]. 

The organization of this paper is as follows. In the next section, we give the main formulas 
of the IQMDD model under the mean field approximation. The main formulas of neutron 
stars in also included. In the third section, some numerical results are presented. The last 
section contains a summary and discussions. 
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II. FORMULAS OF THE IQMDD MODEL WITH p MESON 



A. IQMDD model for nuclear matter 

The Lagrangian density of IQMDD model with a,uj, and p mesons is : 

-U{a) - \f,,F^^ + \mluj,uj^ + Imjp; • pT^ - (3) 

where 

U{a) = \mW + h>a' + \ca' + B, (4) 
-B^ ^al + + '-at (5) 

a, = ^[l + V^-^Kc/b% (6) 

and the quark mass mq{q = u, d) is given by Eq. (1). and are the masses of a and cu 
mesons, F^^ = d^LU^ - d^uj^, = d^'])^ - (9^7?^, g^, gl and gf^ are the couphng constants 
between quark and a meson, quark and uo meson and quark and p meson respectively. 
The equation of motion for quark field under MFA in the whole space is 

[i^-d- (m, - gla) - 7°(^> + ^^^.p)]V^ = (7) 

where a, uj are the mean field values of the a field and the corresponding time component of 
uj field respectively, p is the mean field value of the time component in the third direction of 
isospin for p field, is the third component of the Pauli matrix. The effective quark mass 
m* is given by: 

ml^mq- gla (8) 

In nuclear matter, three quarks constitute a Freidberg-Lee soliton bag [15], and the effective 
nucleon mass is obtained from the bag energy and reads: 

4 ^ I dN 

M*^ = E,E, = ^mf + k\^)dk (9) 
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where 7^ is the quark degeneracy, Kp is Fermi energy of quarks. dNg/dk is the density of 
states for various quarks in a spherical cavity. The expression of dNg/dk adopted in this 
paper con be found in Ref. [16]. 

The Fermi energy Kp of quarks is given by 

3 = ^ttR^ub (10) 



where ub satisfies 



^B^^.V^J C-^)dk (11) 



(27r)3 7o ' dk 

The bag radius R is determined by the equihbrium condition for the nucleon bag: 



5M% 



(12) 



5R 

In nuclear matter, the total energy density and pressure density read 

4„J,.^lte3^1,,.,_|^. (13) 



and 

Pmatter = ^TTT^i I + I ) / ■dp'' + TTTPb 

3 (27r)3 Jo Jo Jm*/+p^ 

2 3 4 8m2^^ ^ ^ 

where 7jv = 2 is degeneracy of proton or neutron, Kp and Kp is Fermi momenta of proton 
and neutron, and pa is the difference between the proton and neutron densities, respectively. 
Therefore 

Pp=^2K'F\Pn = ^,Kf, (15) 

Where Pp and p„ is the density of proton and neutron respectively, and the density of nuclear 
matter ps resds 

PB= Pp + Pn (16) 
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In Eqs. (13, 14), and gp are the coupling constants between nucleon and cu meson, 
nucleon and p meson respectively. They satisfies gu> — and gp — gf, [13]. In MFA, the u 
is determied by baryon number conservation 

QujpB 



(17) 



As that of the QMC model [13], a and p are given by the thermodynamics conditions: 

' d^matter 



/ '^'^ matter \ 



0, and 



dp 



-)r, 



Pb 







respectively. Therefore, p is expressed by 
and a is given by 

2- , ,_2 , -3 

m^a + ba + ca = -— — 

{2ny 



9p 
2ml 



P3 



+ 



N 



-d pi^U 



(18) 



(19) 



(20) 



Eqs. (13-20) form a complete set of equations and we can solve numerically. Our numerical 
results will be shown in the next section. 

Noting that the left hand side of Eq. (20) is a cubic order function of a , except one 
solution a = 0, there are still two solutions. This is a general character for adding a nonlinear 
scalar a field and using MFA to a physical model [17]. But as was pointed in Ref. [1], we 
can prove that one of these solutions corresponds to unstable and unphysical branch, and 
the other corresponds a stable soliton solution. Hereafter we give up the unphysical solution 
and consider the physical solution only. 

We note that the expression for the total energy density, Eq. (13), is very similar to that 
of QHD-II model and QMC model. The differences comes from the effective nucleon mass 
Eqs. (9) and (1), and the self-consistency condition for a field, Eq. (20). Let us consider 



the self-consistency condition and 



[13], we find that the 



,dM 



jR further. Using the same argument as that of Ref. 
g- )ii can be expressed as 

/ 1 \ / 



^ da ' 



R 



-9a X 



1 

Coia) 



for 



I 



QHD - II 
QMC 
y IQMDD J 



model 



(21) 



where the expression of scalar density factor Ci{a) for QMC model can be found in Ref. 
[13]. For IQMDD model, C2(a-) can be obtained numerically. The curves of Ci{a) and C2(a-) 
will be shown in Sec. 3. 
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B. IQMDD model for neutron star 



We now turn to investigate the neutron matter and neutron star for the IQMDD model. 
Since the aim of this paper is to compare the IQMDD model and the QMC model, we use the 
same approximation to study the neutron star as that of the QMC model [14]. More detail 
treatment of neutron star such as the phase transition for the quark matter and the neutron 
matter, the contribution of hyperon and etc, has been neglected. Two basic assumptions: 
the neutron star matter is charge neutrality and reaches to the /^-equilibrium, are adopted 
[14]. Since we assume that the nucleons and light leptons exist in the neutron star only, 
charge neutrality is expressed as 

= E (22) 

l=e,iJ, 

where pi is the number density of particle i{— p, e, /i). Under /3-equilibrium, the processes 

n — >■ p + e~ + Pe, P + e~ n + Ue (23) 

occur at the same rate. This condition can be satisfied when the chemical potentials before 
and after the decay are same. The chemical potential of each particle reads 



/x„ = ^jKf + + _ ]^g^p (24) 



Hp ^ ^JKI + my + g^u + -QpP (25) 



lii^^Kf + mf (26) 

where Ki is the Fermi momentum of the lepton l{e,iJ,). The chemical equilibrium condition 
is expressed as 

/^n = IJ'p + IJ'e, (27) 

/^e = l^ix (28) 

Once the solution has been found, the equation of state(EoS) can be calculated from 

+ 1 5^=^ + \ca^ + ^P3^ + i 5: p^^l^'dK (29) 
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-S^'-k-sif^^-siEf (30) 
Using the Oppenheimer and Volkoff equation 

dM{r) = A7rr^e(r)dr (32) 

where G is gravitational constant and C is the velocity of light, and the equation of state 
for neutron matter given by Eqs. (29), (30) and (9), we can study the physical behavior of 
neutron star for IQMDD model. 

III. NUMERICAL RESULT 

Before numerical calculation, let us discuss the parameters in IQMDD model. First, we 
choose = 783 MeV, rup = 770 MeV and = 509 MeV as that of Ref. [18]. Fixing 
the nucleon mass — 939 MeV, we get B — 174 MeV fm~^. Obviously, the behaviors 
at the saturation point must be explained for a successful model. It reveals that nuclear 
matter saturates at a density po = 0.15 fm~^ with a binding energy per particle E/A = —15 
MeV at zero temperature, and the compression constant to be about K{po) = 210 MeV. 
Therefore we fixed — 2.44, — 4.67, b — —1460 MeV to explain above data. In addition, 
the symmetry energy coefficient asym satisfies 

_ 1 d^e/p) _,9p.2kl kl 

O-sym — 7TV — — )ot=Q — \ ) , r, o / = W-^J 

2 da^ rUp 127r2 G^f^^Mi' 

where 

a=^^^^,A;o = 1.42/m-^ (34) 

Pn + Pp 

Using the data asys=33.2 MeV we fix gp — 9.07. The parameters used to calculate and 
the results of K and for IQMDD model are shown in Table 1. For comparison, we also 
show the corresponding parameters for QMC model in Table 1. The data and results for 



8 



QMC model are taken from Ref. [13]. We see in Table 1 all parameters and results are very 
similar for these two models. Their differences are not remarkable. 

TABLE 1. Comparison of calculated quantities in IQMDD and QMC model 

gl 9l 9l R(fm) K(MeV) M*M{Me V) 

QMC 5.53 1.26 8.44 0.80 200 851 

IQMDD 4.67 2.44 9.07 0.85 210 775 

Our results for symmetric nuclear matter and neutron matter are shown in Fig. 1-3. The 
scalar density factor C{a) as a function of a is shown in Fig. 1 where the dashed curve 
refers to Ci(o") and solid curve to C2(^) respectively. This factor plays an essential role to 
demonstrate the main character of quark structure for different models. We see from Fig. 1 
that the scalar density factors Ci{a) and C2(a-) are both smaller than unity (QHD-II model) 
and decrease when a increases. In particular, C2(o-) is located between the line of unity and 
the curve of Ci(o"). It means that the values of main physical quantities given by IQMDD 
model will almost located between the values given by QHD-II model and by QMC model. 
Our results confirm this conclusion. 

In Fig. 2, we draw the curves of energy per baryon vs. baryon number density for 
both symmetric nuclear matter and for neutron matter respectively. We see that ignoring p 
meson coupling yields a smaller bound state around ps ~ 0.10/m~^ ~ 0.66po in the dotted 
curve, but it becomes unbound solid line when the p meson contribution is introduced. The 
saturation curve for symmetric nuclear matter is shown by dash-dotted curve in Fig. 2. 
In fact, the behavior of curves in Fig. 2 is very similar to that of QMC model, but the 
corresponding value of pb is 0.60po for QMC model and 0.66po for IQMDD model. (See Fig. 
2 of Ref. [13]) 

The equation of state for neutron matter is shown in Fig. 3 where the dashed curve 
presents the result when the p meson contribution is ignored and the solid curve corresponds 
to the full calculation. We see the contribution of p meson is important for the EoS. After 
comparing with the results of QHD-II model and QMC model, we come to a conclusion that 
the shape of equation of state for neutron matter in IQMDD model is qualitatively similar 
to that of QHD-II model and QMC model, it is softer than that of QHD-II model but harder 
than that of the QMC model, as is indicated in Fig. 1 by the behavior of scalar density 
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factor. 

Having shown the IQMDD model can provide an successful description for nuclear and 
neutron matter, we would like to study the structure and composition of neutron stars for 
IQMDD model. Wc will show that it can successfully describe the neutron star. 

The EoS is given by Eqs. (9), (29) and (30) for IQMDD model when the neutron star 
matter reaches to (3 equilibrium. The curve of EoS is shown in Fig. 4. In Fig. 5, we show the 
particle population including e, for different density by solid(n), dashed(p), dotted(e), 
dash-dotted (//) curves respectively. The mass of neutron star in units of sun mass M/Mq 
as a function of central density Sc is plotted in Fig. 6. In Fig. 7 we show the mass radius 
relation of the neutron star. The maximum mass of neutron star M„iax found in Fig. 7 is 
1.73 Mq. It is smaller than the value of 2.2Mq given by QMC model [14]. We would like 
to emphasize that the above treatment for neutron star is too rough. Therefore, the value 
of M^ax is not important. Our aim is to demonstrate that all curves shown in Fig. 4-7 are 
in agreement with those given by QMC model qualitatively [14]. We come to a conclusion 
that perhaps the IQMDD model is a good candidate to replace the QMC model. 



IV. SUMMARY AND DISCUSSION 



In summary, we have added the p meson to the IQMDD model to study the asymmetric 
nuclear matter, especially, the neutron matter and neutron star. The u, d, quarks, nonlin- 
ear scalar a meson field, cu meson field, p meson field and the corresponding quark mesons 
couphngs are including in the new IQMDD model. The isospin effect has been considered 
by introducing isovector p mesons in this model. After fixing the parameters by the ex- 
perimental values such as the massed of nucleon, a meson, uj meson, p meson; saturation 
point, compression constant and the symmetry energy, under MFA, we have investigate the 
physical properties of nuclear matter and neutron matter. We found that the results given 
by IQMDD model are similar to that of QMC model. The values of the main physical quan- 
tities for neutron matter and nuclear matter given by IQMDD model locate in the regions 
between the values given by the QHD-II model and by the QMC model. Employing the 
IQMDD model, we have studied the neutron star and also found its properties almost agree 
with these given by QMC model. We conclude that the new IQMDD model with p meson 
is successful for describing the nuclear matter and neutron matter. Perhaps it can play the 
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role to replace the QMC model. 
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FIG. 1: Scalar density factors C2{cr) and Ci{a) as a function of a for IQMDD model and QMC 
model respectively. 
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FIG. 2: Energy per nucleon for symmetric nuclear matter and for neutron matter. The dash-dotted 
curve is the saturation curve for nuclear matter. The solid curve(with p meson) and the dotted 
curve(omit p meson) show the results for neutron matter. 
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FIG. 3: The equation of state for neutron matter. The dashed hue show the result for omitting p 
meson, while the solid line corresponds to full consideration. 
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FIG. 4: Equation of state for (3- equilibrium neutron star matter. 
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FIG. 5: Populations in neutron star matter as a function of density. 
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FIG. 6: Neutron star mass as a function of the central density. 
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FIG. 7: Masses of neutron stars vs. their radii. 
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